A quantitative study of spin-flip co-tunneling transport in a quantum dot 
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We report detailed transport measurements in a quantum dot in a spin-flip cotunneling regime 
and quantitatively compare the data to microscopic theory. The quantum dot is fabricated by 
lateral gating of a GaAs/AlGaAs heterostructure, and the conductance is measured in the presence 
of an in-plane Zeeman field. We focus on the ratio of the nonlinear conductance values at bias 
voltages exceeding the Zeeman threshold, a regime that permits a spin fiip on the dot, to those 
below the Zeeman threshold, when the spin flip on the dot is energetically forbidden. The data 
obtained in three different odd-occupation dot states show good quantitative agreement with the 
theory with no adjustable parameters. We also compare the theoretical results to the predictions 
of a phenomenological form used previously for the analysis of non-linear cotunneling conductance, 
specifically in the determination of the heterostructure g factor, and find good agreement between 
the two approaches. The ratio of nonlinear conductance values is found to slightly exceed the 
theoretically anticipated value and to be nearly independent of dot-lead tunneling coefficient and 
dot energy level. 
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Electronic transport in nanoscale devices^"— has been 
of significant recent interest, in part for its use as a 
spectroscopic tool for precision studies of fundamental 
phenomena, and because of the relevance of these de- 
vices to spintronics and quantum computationi^"— For 
spintronics, it is important to understand how the spin 
state of a nanosystem couples to its host surroundings. 
Spin-dependent transport can be conveniently studied in 
tunable quantum dots (QD)s»ii"— Using a dot weakly 
coupled to the "leads" with an applied in-plane mag- 
netic field, Kogan et alJ^ showed that the differential 
conductance G = dl/dVds exhibits steps at Vds values 
given by the ratio of the Zeeman energy and the electron 
charge, e, and used a phenomenological fit to the trans- 
port data to measure the heterostructure g factorJ^ii^ 
Later, Lehmann and Loss^i developed a microscopic the- 
ory to calculate the conductance through a QD in this 
regime, which included phonon-assisted spin-flip mecha- 
nisms. In this paper, we present extensive transport data 
of a quantum dot in the spin-fiip cotunneling regime and 
compare the results to microscopic theory Importantly, 
we measure all dot parameters needed for the calculation 
of the conductance, which enables a direct comparison 
between the data and the microscopic theory without 
any adjustable parameters, and find excellent quantita- 
tive agreement between the data and theory. 

We present data obtained for three different choices of 
the dot potential defined by the voltages on the confining 
gates, which correspond to three different occupancies of 
the dot. We focus on the ratio of the device conductance 
above and below the Zeeman threshold as a function of 
the tunneling rate and the dot energy. Since the orbital 



part of the wave function of the two Zeeman spin states 
is the same, the tunneling probabilities for each electron 
crossing the dot depend only on its spin and the spin 
of the dot. Therefore, a useful insight can be obtained 
from the ratio of the device conductance above and be- 
low the Zeeman threshold (i.e., when the bias across the 
dot matches the ratio of the Zeeman energy and the elec- 
tron charge). If the coupling to the leads is extremely 
weak (i.e., the tunneling rates between the dot and the 
leads are much smaller than the spin relaxation rate on 
the dot) one might expect this ratio to be approximately 
2: at large biases, there are two possible dot states (the 
ground spin state and the excited spin state) available 
upon the completion of each tunneling event, whereas at 
low biases, the dot has to remain in the ground spin state. 
In practice, however, the spin relaxation rate due to intra 
dot processes is usually very slow, compared to the tun- 
neling rates in transport experiments between the dot 
and the leadsJ^iiS In that regime, therefore, exchang- 
ing spin with the leads is the dominant mechanism of 
the dot spin relaxation. Predicting the device conduc- 
tance in this regime requires a formalism that includes a 
complete set of rate equations, as we use in this paper 
for a single-orbital, spin 1/2 dotJ^ Our calculations and 
measurements both do reveal a nontrivial value for the 
conductance ratio « 2.4, indicating the important role of 
the current leads in providing spin relaxation in the dot. 

Further, we show that this ratio is independent of the 
dot-lead tunneling rate F over approximately one decade, 
0.02 < F < 0.2 meV, but varies slightly with the dot en- 
ergy, exhibiting a slight minimum in the middle of the 
Coulomb blockade (CB) valley. Finally, we compare the 
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FIG. 1: (Color online) Cotunneling process through a Zee- 
man split orbital occurs when the dot is occupied by an odd 
number of electrons in the Coulomb blockade regime. Tak- 
ing spin up to be the lower-energy state, a spin-down electron 
from the lead can tunnel onto and off the dot resulting in non- 
spin-flip cotunneling as shown in (a). When the bias voltage 
exceeds the Zeeman threshold, |eVtis| > A = \g\fiBB, the 
spin-up electron can also tunnel off the dot resulting in spin- 
flip cotunneling as shown in (b). (c) Micrograph of a device 
nominally identical to the one used in this paper. A quantum 
dot is created after applying negative voltages on electrodes 
Vt, Vs, Vb, and Vq- The electrode Vs is used primarily to 
vary the dot-lead tunneling rate F while the plunger gate Va 
is used to tune the dot energy. Differential conductance is 
measured through source(S) and drain(D) via standard lock- 
in techniques with 2 ^Vrms excitation at 17 Hz. (d) Dif- 
ferential conductance as a function of drain-source voltage 
Vds in the cotunneling regime shows lower conductance (Go) 
at |eVtjs| < A and higher conductance (G+) at |eVds| > A. 
Dashed lines are guides for the average conductance values of 
Go and G+. 



shape of the nonlinear conductance as function of the dot 
bias as obtained from our calculations to the predictions 
of a phenomenological form used in earlier workJ^ For 
a given Zeeman energy, we find excellent agreement be- 
tween the two, which means that either method provides 
a valid choice for using cotunneling transport for g-factor 
measurements. For the device used in this work, using 



both methods, we find the g factor to be 0.2073 ± 0.0013. 

The QD we have studied is created by gating a 
GaAs/AlGaAs heterostructurc. Ti/Au electrodes of our 
single electron transistor (SET) are patterned via e- 
beam and photolithography followed by lift-off. The two- 
dimensional electron gas (2DEG) under the electrodes is 
statically depleted to form an electron droplet (i.e., a QD) 
connected on both sides to the electron reservoirs: source 
and drain [Fig. (He)]. We estimate the diameter of the 
QD to be ^ 0.13 fim, which contains tens of electrons. 
From magnctotransport data wc find that the 2DEG has 
a mobility of 5 x 10^ cm^/(Vs) and an electron density of 
4.8xl0^^cm~^ at 4.2K. The device is oriented parallel to 
the magnetic field within ±1 degree, and is cooled in a 
Leiden Cryogenics dilution refrigerator to a base electron 
temperature Tcicc ^ 55mK. We use standard lock-in tech- 
niques to measure the differential conductance through 
the QD. 

Figure [TJd) shows the differential conductance steps 
at source-drain voltages equal to the Zeeman energy of 
the dot. The tunneling between the dot and the leads is 
relatively weak, so that the Kondo effect in this regime 
is suppressed by thermal fiuctuations. In the Coulomb 
blockade (CB) regime, when the QD has an unpaired 
electron in the dot energy level, the spin degeneracy is 
removed by the Zeeman field, and the level splits into 
spin-up and spin-down states. We label the conductances 
below and above the Zeeman threshold as Go and G+, 
respectively. Figures [TJa) and[IJb) illustrate the possible 
tunneling processes: In an elastic event [Fig. [ija)], the 
dot is left in the ground state and the electron does not 
change its energy as it crosses the dot. If the dot is left 
in an excited state [Fig. (Ub)], the electron energy is 
lowered by A. 

To examine the cotunneling conductance and the ratio 
of G-f to Go quantitatively, we arrange three different 
dot configurations: COT I (Vs = -800 ^ -872, Vt = 
-816, Vb = -1151, Vg = -938 ^ -792 mV); COT 
II {Vs = -960 ^ -1025, Vt = -750, Vb = -1090, 
Vg = -795 -> -671 mV); and COT III {Vs = -800 ^ 
-917, Vt = -750, Vb = -1090, Vg = -1246 ^ -1008 
mV) . For each configuration, the dot contains a different 
number of electrons. To tune the dot-lead coupling, F, 
we use a previously developed computer control of the 
dot gate voltagesi^ and adjust the voltages Vs and Vg 
so as to maintain the occupancy of the dot and keep the 
dot energy in the middle of the Coulomb valley. To tune 
the dot energy \Ei — fJ.\, we vary the plunger gate voltage 
Vg while keeping voltages on other electrodes unchanged. 
We focus on the changes of the conductance as well as 
the ratio G+/Go as either the tunneling rate or the dot 
energy is varied. The experiment is performed for all 
three device configurations described above. 

We measure the tunneling rate F = F^ + F^, where 
^L{R) is the tunneling rate from the left(right) lead, by 
examining the shape of the charging peak as we vary 
the voltages on the gates. Figure dja) shows clearly the 
evolution of the Coulomb charging peak width as Vs is 
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FIG. 2: (Color online) (a) Normalized Coulomb blockade 
peaks taken at Vs= -800, -830, and -872 mV (black, red, 
and blue lines respectively) for dot configuration I (COT I) 
clearly show the variation in peak width, (b) Tunneling rate 
r as a function of Vs- Tunneling rate increases as the side 
gate voltage becomes less negative. Inset: Fitting a Coulomb 
blockade peak (dots) to a thermally broadened lorentzian (red 
solid line) gives the corresponding P. Configurations II and 
III show similar behavior (not shown). 



TABLE I: Capacitance ratio ao, used as the energy lever arm, 
extracted from Coulomb blockade diamonds for different Vs 
values. Tunneling rate F and charging energy U shown are 
parameters for the COT III dot configuration. 



Vs (mV) 


ac 


F (meV) 


U (meV) 


-800 


0.027 


0.19 


2.77 


-850 


0.03 


0.06 


2.89 


-900 


0.036 


0.03 


3.11 



varied. To determine F, we fit the CB conductance line 
shape to a thermally broadened lorentzian (TBL) ^^'^^ 
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In this equation, Vq is the gate voltage that corresponds 
to the CB peak maximum, F is the associated tunneling 
rate, aa is the energy lever arm of the dot, and A is a fit- 
ting parameter which is related to the dot asymmetry^^^ 
S = Tl/Tr. The dot asymmetry for each Vs voltage 
setting is obtained from the height of the CB peak^i^^ 
S varies from 4 to 51 in our measurements. A slight de- 
viation of ac due to a possible shifting of the position 
of the dot has been observed and taken into considera- 
tion. Table I lists aa and other dot parameters for three 
different choices of Vs. To assign the corresponding F 
for each V5, we use the average of the tunneling rates 
extracted from the two adjacent CB peaks in the same 
valley F = {TLp+Tfip)/2, where Fip(flp) corresponds to 
the left(right) CB peak. An approximately linear depen- 
dence of F on Vs, and the TBL fitting to a CB peak are 
shown in Fig. El^b). The overall conductance decreases 
with more negative Vs values, as expected, because of 
the reduction in the transmission of the barriers. 
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FIG. 3: (Color online) (a) Differential conductance as func- 
tion of Vds and Vg at B = 9 T for COT III. The cotunneling 
trace of any given F is taken in the middle of the Coulomb 
valley. Dashed lines with arrows indicate the conductance 
threshold across the odd-occupied valley. Dot-dashed lines 
mark the Coulomb blockade diamond edges, (b) Represen- 
tative differential conductance traces taken in the middle of 
Coulomb valley for different Fs of COT II: from 0.072 meV 
(top) to 0.032 meV (bottom), (c) Comparison of the micro- 
scopic calculation with no adjustable parameters to predic- 
tions of a phenomenological form (Eq. 1 in Ref. [l^) used for 
the data analysis shows good quantitative agreement. The 
gap width is twice the Zeeman energy. 



Figure [3^a) shows the characteristic features of the co- 
tunneling conductance in the presence of the magnetic 
field, for a typical valley with an odd-number electron 
occupation. At each gate voltage, a threshold step is ob- 
served; the separation between the steps at positive and 
negative bias is controlled by the Zeeman energy, and 
it is thus independent of the gate voltage. Figure El^b) 
shows representative traces at several different Fs while 
the dot energy is kept in the mid-point of the Coulomb 
valley as described above. 

In order to make direct comparison between theory 
and experiment, we consider a model where transport 
occurs across a quantum dot contacted to two leads, in 
the presence of a spin-flip mechanism due to the coupling 
of the quantum dot to a phonon bath. The Hamiltonian 
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of the system is described bj*iii2£iS 

H = Ho + Htun + Hsp (2) 

where stands for the Hamiltonian of the isolated dot, 
the ideal leads, and the free phonons, 

i?o = ^ Sana + U nfui + ^ eik-nika + ^ liujqnq; (3) 

cr Ikcr q 

here, Ua {nika) is the number operator of the electron in 
the dot (leads) with spin a and U is Coulomb interaction 
between two electrons in the dot with opposite spins. The 
last term in Eq. [3] describes the free phonon bath with 
occupation numbers Ug and energy hujq. 

The hybridization between the dot and the leads is 
described by the tunneling Hamiltonian 

Htun = ^Ik^kada + H.C, (4) 

Ika 

where c\^^ {da) is the fermionic creation (destruction) 
operator of the electron on the leads (dot). Here we have 
assumed that tunnel matrix elements, Vik, are spin inde- 
pendent. Finally, the spin-phonon interaction is modeled 

by 

Hsp = Y{MqxO-x + Mqy(7y){a)_^ + a,), (5) 

where the bosonic operator (og) creates (destroys) 
a phonon in the mode q; Mqx^Mqy are the spin-phonon 
coupling amplitudes and ax, cry are Pauli matrices. 

To calculate the differential conductance G = dl/dV, 
we derive the current which crosses the quantum dot from 
the left (L) to the right (R) lead. The current through 
the dot can be expressed by^^i^^ 

lLR = eY,WLa'.RaPa (6) 

(TO-' 

where e is electron charge, W^a-'Ma is the transition rate 
for an electron tunneling from L lead (spin a') into R lead 
(spin cr), and can in principle take into account elastic, 
inelastic, as well as phonon-assisted elastic cotunneling 
processes; Pa is the occupation number of electrons in 
the dot which is governed by a master equatioiiii 

dPa 

= -JaaPa + 1 era Pa , (7) 

dt 

with the rate jaa- including spin-flip and inelastic co- 
tunneling processes with a current lead, and spin-flip 
processes due to the spin-phonon coupling. In the sta- 
tionary limit, the solution of the master equation is 
given by = [1 + jga/laa]^^ ■ Detailed expressions 
and discussions for the different rates are found in the 
literature 

Direct comparison of the calculated and experimental 
conductance traces, such as those in Fig. U (right panels). 
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FIG. 4: (Color online) Comparison of microscopic calculations 
(dotted lines) and measurements for Go (squares) and G+ 
(circles) of three dot configurations. Left panels: dependence 
on the tunneling rate. Middle panels: dependence on the 
dot energy. Right panels: microscopic calculation (dotted 
line) of differential conductance as a function of Vds shows 
agreement with the experimental measurement (solid line) for 
the conductance near zero bias, but is slightly off at high bias. 



shows their excellent agreement.— Above the threshold, 
the measured conductance exceeds slightly the calculated 
conductance, arising perhaps from a slight bias depen- 
dence in the barrier transmission coefficients and/or the 
increasing important role of other dot levels ignored in 
the model. We point out that spin-phonon interaction is 
expected to reduce the conductance at high bias; more- 
over, the overshoot seen in the data near threshold is 
not expected for the strongly asymmetric quantum dots 
studied here.— Its nature is still unresolved. 

We have also compared the microscopic theory to 
the phenomenological form used by Kogan et alr^ for 
the analysis of non-linear cotunneling conductance. We 
specifically use both approaches to determine the g fac- 
tor of the hetcrostructurc and find excellent agreement 
[Fig. Eljc)] between both approaches. 

Having obtained the dot energy, F, g factor, and the 
dot asymmetry, we now focus on the conductances (Go) 
and (G+) for the three different dot configurations. Fig- 
urc[4](left panels) shows quantitative agreement between 
the predictions and the data, for over two orders of mag- 
nitude in conductance, as F changes. Notice that the 
theoretical curves are not smooth functions of F since 
the asymmetry factor is not the same for each choice of 
F. The ratio of conductances G+/Go « 2.4, however, is 
nearly independent of F for all three configurations [Fig. 

Next, we address the variations of Gq and G+ with 
dot energy. The dot energy is tuned by varying the 
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FIG. 5: (Color online) (a) The G+/Go ratio as function of F. 
Both calculations (lines) and measurements (symbols) show 
that G+/Go ~ 2.4 is nearly independent of the tunneling 
rate F. (b) The conductance ratio as function of dot energy- 
using as reference the midpoint of the valley. Vertical dashed 
lines indicate where charging peaks appear (at half of the 
charging energies) for all three dot configurations. Charging 
energy values for COT I, II, and III are 2.9, 2.0, and 3.0 meV, 
respectively. The ratio reveals a minimum at the midpoint of 
the valley AE = 0, but it slightly increases as the dot energy 
approaches the charging peaks. 



plunger gate voltage Vq, while maintaining \Ei — ^\/T or 
(U — \Ei — /i|)/r > 4, to avoid the dot entering the mix- 



valence regime. We find that the conductance increases 
symmetrically as the dot energy is tuned away from the 
mid-point of the valley [Fig. |4] (middle panels)]. We ex- 
amine the ratio G+ /Go and find that although nearly 
constant at « 2.4, it exhibits a slight minimum at the 
midpoint of the valley; the ratio increases slightly as the 
dot energy approaches the adjacent CB peaks. Figure 
[5fb) again shows good agreement between the calculated 
and measured results. 

In summary, we have presented a systematic study of 
the differential conductance of a quantum dot in the co- 
tunneling regime for three different dot occupancy con- 
figurations. This allowed us to investigate the depen- 
dence of tunneling rate and dot energy on conductance 
and compare the experimental data to microscopic calcu- 
lations. Independent experiments to determine the pa- 
rameters of the dot state were performed so that com- 
parisons could be made without use of adjustable param- 
eters. Wc find overall excellent agreement between the 
calculations of a simple two-spin quantum dot model and 
the measurements. We find that the ratio of the device 
conductance above the Zeeman threshold to that below 
the threshold is nearly independent of the dot-lead tun- 
neling rate and it is only slightly dependent on the dot 
energy, with a value w 2.4, in near agreement with the 
theoretical ratio. The agreement is best in the middle 
of the Coulomb valley and becomes worse closer to the 
charging peaks, possibly due to the role of higher-excited 
states not included in our calculations. 
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